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Dynamically stable periodic rotations of a driven pendulum provide a unique mechanism for 
generating a uniform rotation from bounded excitations. Numerical bifurcation analysis shows that a 
small ellipticity of the driving has two effects. First, the region in the parameter plane of amplitude 
and frequency of excitation where rotations are possible increases with the ellipticity. Second, the 
resonance tongues, which are the most characteristic feature of the classical bifurcation scenario of a 
parametrically excited pendulum, merge into a single region of instability. 
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I. INTRODUCTION 



The driven pendulum is a paradigm example for study- 
ing nonlinear dynamics in mechanics [1] and beyond 
[2, 3, 4, 5, 6]. Its geometric nonlinearity can be modelled 
reliably (in contrast to, for example, nonlinear damping) 
and mechanical pendula are amenable to experimental 
investigations. The dynamical properties of the classical 
parametrically (harmonic up-down) excited pendulum, 
such as resonances, escape from potential well, symmetry- 
breaking, and periodic and chaotic attractors, have been 
explored in detail experimentally [1, 7, 8, 9, 10] and theo- 
retically [11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 
24, 25, 26]. 

This paper studies what happens to the well-studied 
bifurcation scenarios of the parametrically excited pendu- 
lum if the driving of the pivot of the pendulum follows 
a narrow upright ellipse; see figure 1. Our motivation 
for studying elliptic excitation is that only the elliptic 
component of an arbitrarily shaped periodic excitation 
has an effect on a rotating pendulum for large excitation 
frequencies as encountered in mechanical experiments 
[10]; see section II for a brief explanation. Rotations are 
dynamical regimes of the pendulum with a non-zero long- 
time average angular velocity. Stable periodic rotations 
occur naturally over a large range of excitation parame- 
ters in the parametrically excited pendulum [11]. Thus, a 
rotating pendulum provides a unique mechanism for gener- 
ating a uniformly one-directional rotation from a bounded 
excitation. This is a potential physical principle for har- 
nessing the energy of vibrations, which are not necessarily 
purely in the vertical direction. The other motivation 
for focusing on rotating attractors is that the rotating 
pendulum is ideal for developing and testing non-invasive 
bifurcation and chaos control methods [27, 28] in a real 
experiment: periodic rotations are reliably controllable 
by superimposing feedback control onto the excitation 
without changing the shape of the excitation. This is not 
true in general for oscillations around the hanging-down 
state [29]. 

The dimensionless equation of motion for the elliptically 



excited pendulum is 

(9 + 7^ + (1 + pcos (ut)) sin l9 + epsin(cjt) cos l9 = (1) 

where 7 is the coefficient of linear viscous damping, p is 
the scaled excitation amplitude, uj the rescaled excitation 
frequency, and e is the ratio between the horizontal and 
the vertical diameter of the ellipse traced out by the pivot 
during each period. The classical parametrically excited 
pendulum corresponds to the setting e = 0. 

The two main effects of a small non-zero ellipticity e of 
the excitation are: 

1. The classical resonance tongues for the 1:2 and 
the 1:1 resonance of the parametric pendulum [11] 
merge into a single region of instability of the small- 
amplitude period-one libration (which is a pertur- 
bation of the hanging-down state of the parametric 
pendulum). 

2. The ellipticity breaks the reflection symmetry, which 
causes a preference for rotations that have the same 
direction as the motion of the pivot. Effectively, 
the range of possible excitation frequencies and am- 
plitudes where rotations are supported increases 
for increasing ellipticity. The preferred direction of 
rotation has the same orientation as the motion of 
the base around the ellipse because this rotation 
picks up energy from the additional excitation into 
the horizontal direction. 

These two observations are, in short, the key message 
of the paper. Point 2 is, for large frequencies, universal 
for all shapes of excitation that have a dominant vertical 
component. Section II motivates why we study exci- 
tation along upright ellipses and gives an approximate 
expression for the onset of rotations that is valid also for 
general periodic excitations in the parameter region of 
high-frequency forcing. Section III A and section IIIB 
give two-parameter overviews of the changes to the classi- 
cal structure of resonance tongues and to the existence 
regions of rotations. Section IIIC shows one-parameter 
bifurcation diagrams (for increasing forcing amplitude) 
to illustrate how the different attractors are connected by 
the bifurcations shown in the figures 2 and 3. 
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II. MOTIVATION OF THE ELLIPTIC MODEL 

The angle ^ of a periodically driven pendulum follows 
the equation of motion 



O^jO^fyit) sin i9 + /^(t) cos i9 = 



(2) 



where 7 is the dimensionless viscous damping coefficient 
and the two coordinates of the force on the pendulum bob 
caused by the displacement of the pivot, fx{t) and fy{t)^ 
have period 27r/uj. After rescahng and shifting time, and 
using the convention that ^ = points downwards, the 



ffist Fourier components fx^i and 
have the form 



of fx{t) and fy{t) 



fx,i{t) = epcos{ujt - a), fy,i{t) =pcos{u;t) + 1 (3) 

for a certain p > and a G [0, 27r]. If we assume that the 
vertical component of the forcing is dominant then |e| is 
signicantly less than one. For a large forcing frequency 
uj and correspondingly fast rotation 6{t) in the positive 
direction, the first-order averaged equation of (2) for 
= 0{t) - uot is: 



p 



esina] sin^ + ecosacos^] = 0. (4) 



Positively directed periodic rotations of period one corre- 
spond approximately to equilibria of (4). If we replace a 
by —a in (4) the equilibria of (4) correspond to rotations 
in the negative direction. Approximation (4) implies that 
stable period-one rotations are born in a fold bifurcation 
at 



27CC; = p\/ 1 + ^ 2esina. 



(5) 



Hence, only the first Fourier mode of the excitation de- 
termines the onset of stable rotations in the regime of 
large forcing frequencies. A comparison with averaging 
along the integral curves of the unforced and undamped 
pendulum [15] shows that the high frequency approxima- 
tion is already qualitiatively correct for a; > 1 if all other 
parameters in (2) are of order one and fx{t) and fy{t) are 
not strongly non-harmonic. 

The parameter a determines the inclination of the 
ellipses ranging between excitation along a straight line 
{a = 0) with almost no symmetry breaking effect for 
rotations in (5) and the family of upright ellipses {a = 
±7r/2). 

The extreme case of parametric excitation corresponds 
to (e,a) = (0,0). The other extreme case of horizontal 
excitation, which would be singular (e 00) with our 
choice of parameters, has been studied theoretically by 
Refs. [30, 31]. 

The approximate expression (5) for the onset of rota- 
tions implies that the effects of breaking the reflection 
symmetry (0^0) (— 6>, —0) on rotations are most promi- 
nent if the parametric excitation is perturbed into an 
upright ellipse {a = 7r/2 in (3)). Thus, we restrict our 
study to the specific model (1). Figure 1 shows a me- 
chanical representation of this model: a pendulum having 




FIG. 
lum. 



1: Schematic of an elliptically excited mechanical pendu- 



mass m and length / is driven by a slider on an elliptic 
kinematic constraint. The slider is connected to the pen- 
dulum rod via a pin joint. (In fact, elliptic excitation is 
typical if the pendulum base is floating on water waves.) 
The non-dimensional parameters and the non-dimensional 
time of model (1) can be obtained from the corresponding 
quantities of the mechanical representation by the scaling 
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(6) 



In (6) Y is the amphtude of the vertical displacement ex- 
citation —Ycos{Qt)^ X is the amplitude of the horizontal 
displacement excitation — Xcos(Qt), Q is the driving fre- 
quency, g is the acceleration due to gravity, coq is the linear 
natural frequency of the pendulum at the hanging-down 
angle 6> = 0, / is the length of the (mass-less) pendulor 
arm, m is the mass of the pendulum bob, and c is the vis- 
cous damping coefficient in the mechanical representation 
shown in figure 1. 



HI. RESONANCE STRUCTURE 

In this section we analyse how the introduction of 
a nonzero ellipticity e changes the resonance structure, 
constructing two-parameter bifurcation diagrams in the 
(cj,p)-plane. We also use cross-sections of these diagrams 
(one-parameter bifurcation diagrams) at constant frequen- 
cies CO to make the connection between the different dy- 
namical regimes visible. Throughout our study the bi- 
furcation parameters are the excitation frequency cj, the 
excitation amphtude p and the ellipticity e of the excita- 
tion. The parameter e unfolds the reflection symmetry of 
the parametrically excited pendulum in a symmetric fash- 
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ion: if {0{t),0{t)) is a solution for e then (-6>(t), -0{t)) is 
a solution for — e. Thus, for any ellipticity e the bifurca- 
tions obtained for e are identical to those obtained for — e. 
This implies that we can restrict our attention to e > 0. 

We slice the three-dimensional (cc;,p, e)-space along 
three two-dimensional planes by constructing three two pa- 
rameter bifurcation diagrams in the (a;,p)-plane for three 
different values of e: e = 0, e = 0.1 and e = 0.5. The case 
e = is the classical parametrically excited pendulum 
as studied in [11, 14, 20]. The case e = 0.1 shows how 
a small perturbation of the reflection symmetry affects 
the classical bifurcation scenario and e = 0.5 provides a 
picture of how the bifurcation scenario changes as the 
system deviates further from the parametric pendulum 
case towards a circular excitation. 

The only remaining parameter in the non- 
dimensionalized equation (1) is the viscous damping 7. 
We choose 7 = 0.1 to make our results comparable with 
the results of the previous studies [11, 14, 20]. 

In the following we will discuss the bifurcations of 
oscillations and rotations separately. Oscillations are 
periodic orbits that stay in the potential well around 
the hanging-down position ^ = 0. The average of the 
angular velocity over one period of an oscillation is zero. 
Rotations leave the potential well and have a non-zero 
average angular velocity along one period (for period- 
one rotations the average angular velocity is ±0;). We 
will present rotations and oscillations always in separate 
figures because they coexist over large parameter ranges 
and there is no local bifurcation linking the two types of 
periodic orbits. 



A. Overview of oscillations in (a^,p) plane 

Figure 2 shows the two-parameter bifurcation diagrams 
in the (cj,p)-plane for oscillations. Panel (a) shows the 
classical diagram for e = 0, panel (b) presents the diagram 
for e = 0.1, and panel (c) depicts the diagram for e = 

0. 5. The symbols between bifurcation curves in figure 2 
indicate which attractors are observable in the different 
regions. The most prominent feature of the classical 
diagram 2(a) are the two main resonance tongues where 
the hanging-down position = loses its stability: the 
1:2 resonance at a; ~ 2 and the 1:1 resonance at a; ~ 

1. The 1:2 resonance tongue is bounded by the blue 
period doubling curve H, and the 1:1 resonance tongue, 
which starts at a larger value of forcing {p ^ 0.7), is 
bounded by a pitchfork bifurcation curve (light green 
curves K12 in figure 2(a)). Both tongues are separated 
by a region of stability of 6 = between the curves K2 
and H in figure 2(a). The period doubling H bounding 
the 1:2 resonance has a degeneracy at the point C: it is 
supercritical to the right of C and subcritical to the left 
of C. 

The most significant change for nonzero ellipticity e is 
that the two resonance tongues merge into a single region 
of instability. The period doubling curve H merges with 
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FIG. 2: Bifurcation diagrams in the (a;,p)-plane for (a) e = 0, 
(b) e = 0.1, (c) e = 0.5, computed with (R)AUTO [32, 33]. 
The figure 4 shows the one-parameter bifurcation diagrams 
along the dashed lines in (a) and (b) for cj = 0.87 and a; = 1.8. 
The damping 7 is 0.1. 



one of the non-symmetric period doubling curves L of the 
1:1 tongue. This period doubling (still called H in figures 
2(b) and 2(c)) and the fold curve Ki form the stability 
boundary for the small-amplitude hbration of period one 
around = 0, which is a perturbation of the hanging- 
down position ^ = of order e. The period doubling is 
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The values = 1.8 and uj = 0.87 for these parameter 
paths are representative for the 1:2 and the 1:1 resonance, 
respectively. They are the same as in [20], which studied 
the parametric case e = 0. 



B. Overview of rotations in {co,p) plane 

Figure 3 shows the bifurcations of period-one rotations 
for e = 0.1 (figure 3(a)) and e = 0.5 (figure 3(b)). The 
bifurcations of rotations in the parametric case e = are 
included (J and G, in black and grey) in both panels to 
show the effect of the nonzero ellipticity e. In all cases the 
stable period-one rotations are born, for increasing forcing 

in a fold bifurcation (curves J, and in figure 3) 
and lose their stability in a period doubling (curves G, 
G^ and G^ in figure 3) as p increases further. For even 
higher forcing the period-one rotation regains its stability 
(in the period doublings E, E^ and E^). 

The most notable effect of the nonzero ellipticity is 
that all bifurcations are shifted toward lower forcing for 
rotations in the negative direction (that is, in the same 
direction as the motion of the base along the ellipse). 
The bifurcations of rotations in the positive direction are 
shifted upward. Approximation (5) estimates this effect 
in the limit of high frequency. 

According to figure 1 the rotation in the negative direc- 
tion (0 < 0) rotates in the same direction as the base of 
the pendulum, corresponding to a = — 7r/2 in (5). Thus, 
for e > the curve is shifted downwards from J by 
2juje/{l + e) and the curve is shifted upwards from J 
by 2juje/{l — e) in the limit Cc; ^ oo. The equilibria of 
the averaged equation (4) show that negative rotations 
(a = — 7r/2) pick up energy from the horizontal compo- 
nent of the forcing on average (that is, the factor in front 
of sin (j) is larger than one) whereas the positive rotations 
lose energy. 



One- parameter diagrams for varying forcing 
amplitude 



FIG. 3: Bifurcation diagrams in the (a;,p)-plane comparing the 
bifurcations of the period-one rotations for e = 0.1 (panel(a)) 
and e = 0.5 (panel (b)) with the classical case e = 0, computed 
with (R)AUTO [32, 33]. The figure 4 shows the one-parameter 
bifurcation diagrams along the dashed lines in (a) for uj — 0.87 
and UJ — 1.8. The damping 7 is 0.1. 



subcritical between the points Ci and C2 along the curve 
H. 

The one-parameter bifurcation diagrams along the pa- 
rameter paths marked as dashed lines in figure 2(a) and 
(b) are discussed in detail in section IIIC. They show 
how the other bifurcation curves in figure 2 form the 
stability boundaries for the more complex oscillations. 



Figure 4 shows a series of four one-parameter bifurca- 
tion diagrams for varying forcing amplitude p. We pick 
two values for the frequency (the same as in [20]): a; = 1.8 
(panel (a) and (b)), which is in the 1:2 resonance tongue, 
and uo = 0.87 (panel (c) and (d)), which is in the 1:1 
resonance tongue, and two values for the ellipticity: e = 
(panel (a) and (c)) and e = 0.1 (panel (b) and (d)). The 
vertical axis of all panels shows the coordinate of the 
stroboscopic map of (1). Stable oscillations and rotations 
are dark green thick lines, unstable oscillations and ro- 
tations are bright red thin lines. All bifurcation curves 
in the two-parameter diagrams figure 2 and figure 3 have 
been constructed by continuing the bifurcations shown as 
dark circles in figure 4. We have shifted the value of by 
27r for all rotations to avoid the curves illustrating bifurca- 
tions associated with rotations and oscillations obscuring 
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FIG. 4: One parameter bifurcation diagrams for the strobo- 
scopic map ^ for different values of u;, varying the scaled 
forcing amplitude p as the bifurcation parameter. The other 
parameters are e = 0, a; = 1.8 in panel (a), e = 0.1, cj = 1.8 in 
panel (b), e = 0, a; = 0.87 in panel (c), and e = 0.1, cj = 0.87 
in panel (d). Black dots are the attractors computed for the 
initial condition Oq = O.OItt, = with Dynamics [34]. The 
stable (thick dark green) and unstable (thin light red) periodic 
orbits, computed with (R)AUTO [32, 33], are superimposed. 
Oscillations are shown in the ^-range of [— 7r,7r], rotations are 
shown in the ^-range [7r,37r]. The damping 7 is 0.1. The 
basins of attraction shown in figure 5 have been computed at 
p = 0.5, p=l,p=1.4 (see vertical dashed lines in panel (a) 
and (b)). 



each other. The underlying black dots show the long-time 
behavior of the initial conditions (6>o, ^0) = (O.OItt, 0) (and 
(^0, ^0) = (27r + O.OItt, 0)) after waiting for a transient of 
1000 periods of excitation, computed with Dynamics [34]. 



figures 4(a) and (c). The same apphes to the bifurcations 
of the nonsymmetric orbits: the fold J and the period 
doubling G of the rotations, and the period doublings F 
and L for the nonsymmetric oscillations (starting rapidly 
accumulating period doubling cascades) are symmetric 
pairs of bifurcations, occuring simultaneously. This sym- 
metry is broken by the increase of e such that formerly 
symmetric branches are now different: rotations in the 
negative direction emerging from the fold J" already exist 
for smaller forcing p than the rotations in the positive 
direction emerging from Jp, which are shifted toward 
larger forcing p. Similarly, the formerly symmetric pairs 
of nonsymmetric oscillations lose their symmetry: one 
family is always shifted toward larger p (born at the fold 
F2 in Figure 4(b), and K2 in Figure 4(d)), the other 
family becomes a continuous extension of the formerly 
symmetric oscillation. 

For a; = 1.8 the visibility of chaotic attractors (bands of 
small black dots in figure 4) is shifted toward larger p by 
the symmetry breaking because stable periodic rotations 
exist for larger p (up to p ~ 1.35). At a; = 0.87 the 
simulation also showed period-two oscillations jumping 
between two potential wells for larger p in the simulation 
results (black lines evident after the chaotic bands in 
panel (c) and (d)). 

Figure 5 shows how the basins of attraction lose their 
symmetry when one increases e from to 0.1. The color 
coding of each point in the (^, ^)-plane is chosen according 
to the attractor which the stroboscopic map reaches start- 
ing from this point. The forcing is p = 0.5 in figure 5 (a) 
and (d), p = 1 in panel (b) and (e), and p = 1.4 in panel 
(c) and (f). The periodic attractors (shown as dots) with 
^-coordinates ^ ±2 correspond to periodic rotations. The 
panels (d) and (e) show that for small p {p = 0.5 and 
p = 1) one direction of rotation (negative) has a visibly 
larger basin of attraction than the other. At p = 1.4 
the change of e caused a crisis of the chaotic attractor in 
panel (c), creating a period-three rotation. 

For further increase of e the effect that one attractor 
of the formerly symmetric pair of nonsymmetric periodic 
orbits is shifted toward higher values of forcing becomes 
more pronounced (as shown in the the two-parameter 
diagrams in figure 2(c) and figure 3(b)). This shift de- 
pends strongly nonlinear ly on e: for example, formula 
(5) already underestimates this shift for Jp by 20% for 
e = 0.5, which still corresponds to a narrow ellipse. 



The main feature of the transition from e = to nonzero 
e is the breaking of the reflection symmetry. The sym- 
metric system (with e = 0) has the pitchfork bifurcations 
A (for period two in Figure 4(a)) and Ki (for period one 
in Figure 4(c)) linking families of symmetric and non- 
symmetric oscillations. These pitchfork bifurcations are 
perturbed into fold bifurcations (also named A and Ki in 
the figures 4(b) and (d)). 

The rotations (which are nonsymmetric orbits) and the 
nonsymmetric oscillations come in pairs of orbits symmet- 
ric to each other and lying on top of each other for e = in 



IV. CONCLUSIONS 

Introducing a horizontal component into the excita- 
tion of the classical parametrically excited pendulum 
results in a symmetry breaking scenario. The excitation 
changes from a purely vertical motion to a motion on an 
ellipse. The main effects of this ellipticity are: first, the 
well-known 1:2 and 1:1 resonance tongues of the para- 
metric pendulum merge into a single region of instability, 
bounded by a period doubling and a fold (saddle-node) 
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FIG. 5: Basins of attraction for different values of the forcing 
amplitude, p, and ellipticity e, computed with Dynamics [34]. 
Panel (a): e = 0, p = 0.5, (b): e = 0, p = 1, (c): e = 0, 
p = 1.4, (d): e = 0.1, p = 0.5, (e): e = 0.1, p = 1, (f) e = 0, 
p = 1.4. The other parameters are 7 = 0.1 and 00 = 1.8. 
Stable periodic orbits are shown as large dots (their color is 
chosen to give a contrast to its own basin of attraction). 



overall forcing amplitudes. 

Both effects of ellipticity are favorable for rotation: 
the first effect implies that small-amplitude oscillations 
around the hanging-down position, which are attractors 
competing with rotations, lose their stability for smaller 
forcing compared with the e = case. The second ef- 
fect means that the parameter region in the frequency- 
amplitude plane where rotations are supported increases 
with increasing elhpticity. The excitation amphtude neces- 
sary to sustain rotations in both directions also increases 
with increasing e because of the increasing imperfection 
of the symmetry. 

A comparison between the bifurcation scenarios of the 
model and the actual experiment is still outstanding. Di- 
rect bifurcation analysis for experiments is a challenging 
task that may require the development of entirely new 
experimental methods. Apart from the experimental ver- 
ification, other open questions are: the high-frequency 
approximation (5) suggests that for a circular excitation 
(e = 1, a = 7r/4) rotations against the base excitation are 
impossible regardless of the level of forcing and damping. 
This is not true in general for lower frequency and suffi- 
ciently small damping. Thus, we expect that, depending 
on the shape of the excitation, there must be a critical 
damping level below which rotations against the excita- 
tion direction become possible for a suitable range of the 
excitation amplitudes and frequencies. 

The small dissipation restricts the type of bifurcations 
and regimes encountered in the system (for example, torus 
bifurcations are impossible). We expect that even a small 
amount of interaction between the pendulum and the 
base will lead to large regions in the frequency- forcing 
plane where one can observe quasi-periodicity. Escape 
from a potential well tends to lead to indeterminacy as 
introduced in [35] . The precise sequence of heteroclinic 
tangencies leading to escape from the potential well is 
still open even for the parametrically excited pendulum. 



bifurcation of the small amplitude oscillation. Second, 
rotations of the pendulum that have the same direction 
as the base motion pick up energy from the horizontal 
excitation such that they are already present at lower 
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